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The eiconal and ray equations for a stratified moving medium are analyzed to determine 
which rays connect a source located at an altitude z = z s to an observer on the ground 
plane. The characteristics of all available rays are ascertained. Further, the rays which 
determine the shadow boundary are found. Thus, the shadow boundary may be easily 
calculated and any observers for which an eigenray does not exist eliminated from further 
consideration. Therefore, the analysis may be applied to save computer time in two 
ways. First, rays that do not reach the ground may be eliminated before they are traced, 
and observers for which no eigenray exists may be eliminated from consideration. 


I. Introduction 

The geometric acoustic theory 10 provides a very flexible method for the study of sound 
propagation in the earth’s atmosphere. However, in the geometric theory the 
determination of the acoustic field at a fixed observer location requires the solution of a 
two-point boundary value problem for the eigenray. If the medium through which the 
disturbance propagates is assumed to be stationary and homogeneous, this problem has a 
very simple geometric solution — so simple, in fact, that we solve it without realizing that 
a problem has arisen — and very little computation time is required to detennine the 
eigenray and the geometric spreading losses. However, if the sound speed is assumed to 
vary, perhaps with altitude alone as would be fitting for atmospheric propagation studies, 
two problems arise which are absent in the stationary, homogeneous medium problem. 
First, a shadow boundary, which divides the plane on which the observers reside into an 
ensonified region and a zone of silence, or shadow zone, may exist. If such a shadow 
region exists, there are no rays connecting the source to observers in that region. Thus, 
any computer time spent attempting to find such a ray is wasted. Second, some rays that 
would reach the observer plane in the stationary, homogeneous medium case, fail to do so 
when the sound speed varies with altitude, while some rays that would not reach the 
ground in the stationary, homogeneous medium case, may do so when the sound speed 
varies with altitude. Clearly, tracing rays that do not reach the observer plane wastes 
computer time, while not considering rays that do leads to erroneous noise predictions. If 
only the stationary inhomogeneous medium case is of interest, simple well known criteria 
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can be applied to overcome these problems 6 ' 7 . If, however, the effects of a sheared mean 
flow are important, the required criteria are neither so simple, nor as well known. 

In the following, analysis of the eiconal and ray equations will provide a complete 
description of the behavior of any ray before that ray is traced, allowing rays that do not 
reach the ground to be eliminated from consideration. Further, the set of rays that make 
up the shadow boundary is detennined. Thus, the tracing of a relatively limited number 
of rays allows the determination of the shadow boundary and the elimination of observers 
that cannot be reached by any ray. 

II. Notation and Fundamental Equations 


In the following the main interest is determining which rays will carry energy from a 
source at an altitude z s to an observer at z = 0. The source is assumed to be embedded in 
a medium whose sound speed c(z) and mean flow velocity 

U(z) = U(z)x + V(z)y = |uj[xCos(/7) + ySin(/?)] are functions of z alone. Here, and 


in the following, a caret over a letter indicates a unit vector. Thus, for example, the 
notation x indicates a unit vector directed in the positive x direction. Similarly, an arrow 
over a quantity indicates a vector that need not be a vector of unit length. In the 
following, consideration will be limited to obtaining useful information from the eiconal 
equation 4 


2 , 2,2 [1 - (US X + VS )] 2 

+ s y + s z = 1 z 


( 1 ) 


and the ray equations 


Note that the vector 



v 

cs + U 


s x x + s y y + s z z 


s: + s! + s 2 


( 2 ) 


( 3 ) 


is a unit vector that is normal to the wave-front, or surface of constant phase, and that 
equation (2) simply states that the surface of constant phase moves nonnal to itself at the 
local sound speed c while being carried along at the local flow velocity by the mean flow. 
The three-dimensional curve x = x r obtained by integrating equation (2) is the path 
followed by a point on the constant phase surface as the disturbance propagates. It is, 
therefore, a ray. Employing equation (1), equation (3) may be written in the alternate 
form 
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[l - (US, + VS y )] 

For the stratified medium to which we restrict our attention, the quantities s x and s y are 
constant along any given ray 4 , and, given s x and s y , s z is determined algebraically by 
solving equation (1). Hence, the rays are a two parameter (s x ,s y ) family of curves in three 
dimensional space, and the following analysis seeks to limit the range of s x and s y such 
that any pair (s x ,s y ) in the restricted range will provide a ray that connects a source on the 
plane z = z s to an observer on the plane z = 0. 


ITT. The Eiconal Equation as an Equation for a Prolate Spheroid. 


An alternate form of equation (1), 


, MA + M ySy _ 

M ’ 1 


~ M tfy M x s y -M y s , 


1 

2 ' + - 

1 

r + 

1 

L(1 - M 2 )c_ 


_V(1 - M 2 )c 1 


_V(1 - M 2 )c_ 


= 1 (5) 


of value for some of the discussion to follow, is obtained by introduction of the notation 

M = — = M x x + My (6) 

c 

M = -^M 2 + M 2 (7) 

and minor algebraic manipulation. Note that equation (5) is the equation of a prolate 
spheroid. 

Turning Points 

At a turning point the ray is tangent to the plane z = z t . Propagation of the disturbance in 
the z direction is reversed at this point, and it follows from equation (2) that s z = 0. This 
paper essentially reduces to a discussion of the turning points, i.e., those values of s x and 
s y that satisfy the equation 


s; + s; = 


[1 - (U t s x + V t s )] 2 


( 8 ) 
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obtained by setting s z = 0 in equation (1). Here, the subscript t implies that the 
subscripted quantity is evaluated at the turning point. Of particular interest, since they 
form the shadow boundary, are the turning points located at the ground. Since the mean 
flow is assumed to be zero at the ground, equation (8) becomes 




s: + s; = ^ 


(9) 


for this particular single parameter set of turning points. Here and in the following Co is 
the sound speed at the ground. 


IV. A Set of Significant Ellipses. 

Since equation (5) is the equation of a prolate spheroid, the equation 


, m A + M y s y _ 

1 M ' " L 


~ M M x s y -M y s 


1 

2 ' + - 

1 

Led - m 2 )J 


_cV(l - M 2 )_ 


= 1 


( 10 ) 


obtained from it by setting s z = 0, is the equation of an ellipse. That is, all rays that have 
a turning point at a given altitude z possess (s x ,s y ) values that fall on this ellipse, referred 
to here as the ellipse of turning points (ETP). Note that both a coordinate rotation and a 
coordinate translation are required to place the coordinate origin at the center of the 
ellipse and to align the coordinate axis with the principal axis of the ellipse. The rotation 

V 

is through the angle /? = Arctan(— ). Thus, the angle of the major axis of the ellipse to 

the s x axis is equal to the angle of the mean flow vector to the x-axis. The location of the 
center and the foci of the ellipse, and the eccentricity of the ellipse, easily determined 
because equation (10) is in standard fonn, are presented in the following Table. Further, 
it should be noted that for analytic functions M x (z), M y (z), and c(z), the motion of the 
ellipse is continuous as z changes. 


Table. Properties of the Ellipse 


Center of ellipse 

„ _ r -MCos (J3) _ -MSin <J3) 

x c(l - M 2 ) ’" y c(l - M 2 ) 

Foci 

(s x ,s y ) = (0,0); 

- 2MCos(/?) - 2MSin(/?)-| 
1 x ’ y) 1 c( 1 - M 2 ) ’ c( 1 - M 2 ) J 

Eccentricity 

M 
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In order to visualize the variation of the ellipse with altitude, consider an atmosphere 

v 

within which the wind speed, U(z) = U(z)x, increases monotonically with altitude, and 
for which the temperature, hence the sound speed, is constant. Then, the Mach number 
of the atmosphere will increase with altitude, and the set of curves generated either by 
equation (8) or, equivalently, equation (10) will consist of a circle centered at (s x ,s y ) = 
(0,0), obtained for z = 0, and ellipses with a common focus at (s x ,s y ) = (0,0). The second 
focus which, for z > 0 separates from the first, moves in the negative s x direction along 
the s x axis. The distance of this second focus from the first increases with increasing 
Mach number, hence altitude, and, to first order in the Mach number, this distance is 
directly proportional to the Mach number, see the Table, Row 2. Also, see Figure 1, 
which shows three of these ellipses for an atmosphere with a flow speed increasing at 
0.024 (ft/sec)/ft. The first ellipse is for z = 0, the second for z = 2,500 ft., and the third 
for z = 10,000 ft. The motion of the center of the ellipse, and the second focus of the 
ellipse, along the negative s x -axis is clearly indicated. Note that each ellipse has an 
eccentricity equal to the Mach number. Since M « 1, each ellipse is nearly a circle with 
a radius r = 1/c. Thus, if the temperature, hence the sound speed, decreases with altitude, 
as is normally the case, the radius of the approximating circle will increase with altitude. 
If, however, the temperature increases with altitude, the radius of the approximating 
circle will decrease with altitude. 
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sx, sec/ft 

Figure 1. The motion of the second focus as z increases in an atmosphere within which 
the Mach number increases with altitude. The first focus is at (s x ,s y ) = (0,0). 

U = (0.024 z/sec)x, z in feet, c = 1 100 ft/sec. 


V. The Atmosphere as a Ray Filter 

For any value of z, that is for any altitude, equation (10) is the equation of an ellipse in 
the (s x ,s y ) plane. This ellipse divides the (s x ,s y ) plane into two regions. For points (s x ,s y ) 
outside of the ellipse propagation is evanescent in the z direction. The interior of the 
ellipse consists of (s x ,s y ) pairs that define rays that penetrate the plane. The boundary 
separating these regions, the ellipse itself, consists of (s x ,s y ) pairs that define rays that 
turn at this value of z. Thus, one way to look at the atmosphere, at least as to how it 
affects an acoustic disturbance propagating through it, is as a set of these filters, one for 
each plane z = constant. Each filter either removes certain rays from the field at the 
observer, or turns these rays back towards the observation plane. 

If all of the ETPs for a given atmosphere were drawn in the same (s x ,s y ) plane, the plane 
would be divided into four regions, see Figure 2. For (s x ,s y ) in region 1, propagation in z 
is evanescent. These waves either grow or decay as they propagate from the source to the 
receiver. Thus, in the first case, these waves must be removed from the field by a 
judicious application of an initial condition, and they may be neglected in the second. 
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Hence, these values of (s x ,s y ) need not be considered. Any ray specified by (s x ,s y ) in 
regions 2 and 3 will go through at least one turning point. These rays will be considered 
in more detail later. Finally, for any (s x ,s y ) in region 4, the specified ray will connect the 
source to some point on the ground without going through a turning point. Clearly, a 
robust ray tracing program would easily be provided if consideration were restricted to 
rays specified by (s x ,s y ) in region 4. However, in the following, we will take a closer 
look at regions 2 and 3 to show how mapping out these regions allows consideration of 
all relevant rays while easily determining rays that cannot contribute to the field at the 
observer plane. 



sx sec/ft 


Figure 2. The set of all ellipses for a given atmosphere divides the (s x ,s y ) plane into four 
regions. 


VI. The Initial Refinement, Hacking Away at the Boundary 

If, instead of plotting all of the ellipses for a given atmosphere, a minimalist approach is 
taken, and only three ellipses, one each for the observer plane, the source plane, and the 
plane representing the top of the atmosphere, are drawn, a rather skeletal representation 
of Figure 2 is produced. This representation is presented in Figure 3. The three ellipses 
considered here represent the boundaries of significant regions. Since rays associated 
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with (s x ,s y ) in Region 2 of Figure 2 carry evanescent waves, points (s x ,s y ) to the right of 
the ellipse for the source plane may be ignored. This suggests use of the right half of the 
source plane ellipse as part of the boundary of the region of interest. Further, points 
(s x ,s y ) which fall outside the circle for the ground plane provide rays which cannot reach 
the ground. This suggests using the left half of the circle for the ground plane as the 
remainder of the boundary of the region of interest. Taking, then, the intersection of the 

regions s x + s^ < \ and s x + sj < 

L, 0 

plane that more closely approximates the region containing all (s x ,s y ) pairs that determine 
rays that can contribute to the field at the observer plane. This region is shown in Figure 
4, which represents a first refinement of the region suggested earlier. 



Figure 3. A skeletal representation of the elliptical acoustic filters for an atmosphere. 


(U s s x + V s s v )]- 


yields a region in the (s x ,s y ) 
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Figure 4. The region (s x ,s y ) for rays that will certainly reach the ground, and two regions 
requiring further analysis. 


Except for points (s x ,s y ) located within the small areas of concern indicated in Figure 4, 
all points (s x ,s y ) in this region provide rays that will reach the ground. Further, all points 
(s x ,s y ) outside of this region either fall outside of the circle s 2 + s 2 = c^, 2 , and therefore 
determine rays that cannot reach the ground, or fall outside of the ellipse 
£ + 4 = [i - (u.s, + Vs y )] 1 /cl, and, therefore, detennine rays which carry 
evanescent waves at the source. Generally, these may safely be ignored. Flence, 
consideration may be limited to points (s x ,s y ) in the indicated region. 

VII. The Final Refinement 

Now consider the regions indicated as “Areas of concern” in Figure 4. Within these 
regions, it is conceivable that an ellipse associated with an altitude between that of the 
source and the ground passes through our region of interest. The points (s x ,s y ) falling on 
this ellipse would then be associated with rays passing through a turning point. These 
points (s x ,s y ), then, are associated either with ducted rays or with rays that pass through a 
turning point and do not reach the ground. Consider the line segment between the 
origin (s x ,s y ) = (0,0) and an arbitrary point on the border of the region of interest indicated 
in Figure 4. Fet this line segment be at an angle a to the positive s x axis. Then, the 
intersection of any ellipse entering the region of interest with this line segment will be a 
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point. As z varies, 0 < z < z m , where z m is the altitude of the top of the atmosphere, this 
point will move along the line segment, as indicated in Figure 4. The vector from the 
origin to this moving point will be denoted by s t . Thus, the equation 

s t = |s t |[xCos(a) + ySin(a)] (11) 

must hold. Further, since this point lies on an ETP, s z = 0, and it must also be true that 

iV ,2 [1 - |u||s t |Cos(a - p)f 

h = — — 2 ( 12 ) 

c 

There are two values of |s t | that satisfy this equation, but for subsonic flow the only 
physically valid root is 


Is I 


[c + |u|Cos(a - /?)] 


(13) 


Differentiating equation (13) with-respect-to z while holding a constant, and setting 
0, yields 


d|s, 


dz 


dU 

dz 


V 




dz 


(14) 


Elimination of z between equation (13) and equation (14) yields an equation for the 
minimum value of |s t | as a function of a. The curve so generated, shown in Figure 5, is 
the envelope of the ellipses that sweep through the circle s 2 + s 2 = c,', 2 . There is, 

obviously, more to the envelope curve than is shown here. However, all that is really of 
interest are the segments of the envelope curve between where it is tangent to the circle 
S' + s; = and where it is tangent to the ellipse 

s 2 + s 2 ; = 11 - (U s T + V, s )l 2 /c 2 . Here, the subscript z m indicates evaluation of 

x y L v z m x z m J ' z m 

the indicated quantity at the top of the atmosphere. There are two of these segments, one 
each for the two areas of concern indicated in Figure 4. For the case considered here, 
however, these are mirror images of each other. Therefore, only the upper segment is 
shown in Figure 5. 
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Figure 5. The envelope curve obtained by elimination of z between equations (13) and 
(14) connects the circle s 2 + s 2 = c' 0 2 with the ellipse 

s 2 + s 2 ; = 11 - (U s T + V, s„)l 2 c' 2 . Note that each point on this curve is associated 

x y z m x z m y z m 

with a value of z, 0 < z < z m , and that its end points are associated with z = 0, the ground, 
and z = z m , the top of the atmosphere, as indicated in the figure. 

The two points at which the envelope curve is tangent to any ellipse divide that ellipse 
into two segments. Along the envelope curve, hence at the point where the envelope 
curve is tangent to the ellipse, equation (14) is satisfied. Therefore, along one of the two 
segments of the ellipse it must be true that 

dU s t 

~7 ~ ' IV" 
dz |s t 

while along the other the inequality 

dU l 
dz s t 

must hold. Consider a wave approaching the turning point from below, as shown on the 
left-hand-side of Figure 6. In this case, s 2 , while zero at z = z t , must be greater than zero 
in z < z t , at least near z = z t . Generally, in this case, s 2 < 0 in z > z t . Flence, we must 
have d(s 2 )/dz < 0 at z = z t . Differentiating Equation (13) with-respect-to z, and requiring 
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d(s^)/dz < 0, yields inequality (15). Similarly, if we consider a wave approaching the 
turning point from above, as shown on the right-hand-side of Figure 6, we must have 
d(s’)/dz > 0, and it may be shown that inequality (16) must hold. Therefore, the points 
where the envelope curve is tangent to the ellipse divide that ellipse into two segments. 
Along one of these segments the corresponding rays approach the turning point from 
below, with evanescent propagation above the turning point, whereas along the other 
segment the corresponding rays approach the turning point from above, with evanescent 
propagation below the turning point. Figures 7 and 8 summarize these results. 


Waves Approaching The Turning Point 
From Below 


s: < o 


Evanescent Waves 



Waves 


d s z n 

<0 

dz 


dU I dc ~ 

— > 0 


dz s t 


dz 


Waves Approaching The Turning Point 
From Above 

Sz > 0 

Propagating Waves 



ds[ 

dz 

> 0 

z = z t 

dU s t 

dc 

, ■ v 

T + < 0 

dz s t 

dz 


Figure 6. The inequalities (15) and (16) and their implications for the direction of a ray at 
a turning point. 


12 

American Institute of Aeronautics and Astronautics 


44 th AIAA Aerospace Sciences Meeting and Exhibit 
9-14 January 20006, Reno, Nevada 


AIAA-2006-0412 


~TI slahralrw7a. : 

(sx,sy) pairs with turning point at this altitude, 

and on this side of the envelope curve, 

are filtered out at this altitude. For a source below 

this altitude these rays carry an evanescent disturbance 


Turning points 
Envelope curve 
Turning points 



On this segment 
of the ellipse 


dU s, 

dc 


—— -ATT + 


< 0 

dz |s,| 

dz 




(sx,sy) pairs with turning point at this altitude, and on 

this side of the envelope curve, are turned toward the ground at this altitude. 
Therefore, these rays must be ducted, although they may not reach the ground. 
For source above this altitude these rays carry evanescent disturbance. 


-0.0005 


0.0005 


0.001 


0.0015 


Figure 7. How the envelope curve divides any given ellipse of turning points, and the 
implied direction of the ray. 


Ground 


Arc of the Source Plane Ellipse of turning 
points 

Envelope Curve of ETPs 


0.0008835 

0.000143 



Turning Points 
BELOW 
the source 





+ * <0 
dz |s,| dz 




Along this section 
of the ellipse 


FILE: RENO-06-lDB.xls 


Top of the 
atmosphere. 


Figure 8. A closer look at the neighborhood of the point where the envelope curve is 
tangent to the source plane ellipse of turning points. 
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IIX. The Ray Behavior For Each Region 

It is now possible to give a rather complete description of the characteristics of a ray 
determined by any point in the (s x ,s y ) plane. It is suggested that the reader refer to 
Figures 9a-c while reading the following description. Figure 9a provides a close-up view 
of the upper “Region of concern” shown in Figure 4, while Figure 9c provides a close-up 
view of the lower “Region of concern” of that figure. Figure 9b provides a birds-eye 
view of the (s x ,s y ) plane. Note that regions referred to by number are regions seen, in 
some sense, when a rather gross view of the (s x ,s y ) plane is taken. Regions referred to by 
letter are regions that would be seen only upon a rather close look at the (s x ,s y ) plane. 

Referring to Figure 9b, points (s x ,s y ) in Region 1 are outside of all of the ellipses for 
points within the modeled atmosphere. Thus, for these points (s x ,s y ), s z is imaginary for 
all z. Propagation is, therefore, evanescent in z. These values of (s x ,s y ) may be 
neglected. 

For points (s x ,s y ) in Region 2, s z is imaginary at all altitudes between the source and the 
observer. These values of (s x ,s y ) may be neglected. 

For points (s x ,s y ) in Region 3, s z is real at the source, but imaginary at the observer plane. 
Therefore, these rays must pass through a turning point between the source and the 
observer, and do not reach the ground. These values of (s x ,s y ) may be neglected. 

For points (s x ,s y ) in Region 4, s z is real both at the source and at the observer plane. 
Further, except for a small subset of points (s x ,s y ) within this region, and considered 
later, s z is real at all altitudes between the source and the observer. Thus, except for a 
small subset of (s x ,s y ) in this region, these rays carry energy from the source to the 
observer plane. 

For points (s x ,s y ) in Region 5, s z is real at both the source plane and the observer plane, 
and at all points between the source and the observer. Once again, these rays carry 
energy from the source to the observer plane. However, for these points (s x ,s y ), a turning 
point occurs somewhere above the source altitude. Therefore, these rays are ducted, and 
multipathing occurs. 

For points (s x ,s y ) in Region 6, s z , although real at the observer plane, is imaginary at the 
source. These points (s x ,s y ) need not be considered. 

We now take a closer look at the regions of concern of Figure 4. The reader should 
consult Figures 9a and 9c during this discussion. 

Region A is inside both the ground circle and the ellipse for the source plane. Further, no 
turning points have entered this region, as it lies to the left of the envelope curve. Region 
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A is, in fact, Region 4 discussed above. Each ray associated with a point (s x ,s y ) within 
this region carries energy from the source to the plane z = 0. 

Region B is outside of the ground circle, but inside the source ellipse. Therefore, s z is 
real at the source, but imaginary at the ground. Hence, there is at least one turning point 
between the source and the ground plane for points (s x ,s y ) within this region. The ray is 
turned away from the observation plane at this turning point. Therefore, these pairs 
(s x ,s y ) need not be considered. 

Regions C and E are both within both the source ellipse and the circle for the ground 
plane. Thus, s z is real at both the source and the receiver. However, the behavior of the 
ray is different for a ray associated with a point in region C than it is for a ray associated 
with a point in Region E. Referring to Figure 8 it may be seen that the turning point for a 
ray from Region C is below the source, while the turning point for a ray from Region E is 
above the source. Thus, rays from Region C are turned away from the observation plane 
while those from Region E are turned toward the observation plane. Therefore, rays from 
Region C may be ignored while those from Region E are ducted between an altitude at or 
above the source plane and the ground and contribute to the disturbance field. Note that 
for rays associated with points in Region E there are two turning points above the source 
altitude. At the first, the ray is turned back toward the source, hence the ground. At the 
second, a ray is turned away from the ground. 

Region D is within the ground circle, but outside the source plane ellipse. Therefore, for 
points (s x ,s y ) within Region D, s z is imaginary at the source, although real at the observer. 
Again, these pairs (s x ,s y ) may be neglected. 

Region F is within the ellipse for the source plane, and inside of the circle for the ground 
plane. Further, Region F is outside of the envelope curve. Note, however, that Region F 
is also outside of the ellipse for the top of the atmosphere. Rays associated with points 
within Region F go through a turning point at an altitude above the source plane. In 
contrast to points in Region E, however, there is only one turning point above the source 
for rays associated with points in Region F. Region F is, in fact, Region 5. These rays are 
ducted and contribute to the disturbance field on the observer plane. 

For points (s x ,s y ) within Region G, s z is real at the ground, but imaginary at the source. 
These rays may be ignored. 

For points (s x ,s y ) within Region H, s z is imaginary for all points within the atmosphere. 
Thus, Region H is Region 1 . 

For Region I, s z is imaginary for all points between the source and the observer. This is 
Region 2. 

In summary, only rays from Region E, in addition to rays from Regions 4 and 5 
considered previously, will reach the ground and add to the disturbance field. Rays from 
all other regions should be ignored from the outset. 
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Figure 9a. Magnified view of the upper Area of concern from Figure 4 and the regions 
formed in that plane by the ground circle, source plane ellipse, the ellipse for the top of 
the atmosphere, and the envelope curve. 
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Figure 9b. An overview of the (s x ,s y ) plane and the regions formed in that plane by the 
ground circle, source plane ellipse, and the ellipse for the top of the atmosphere. 
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Figure 9c. Magnified view of the lower Area of concern from Figure 4 and the regions 
formed in that plane by the ground circle, source plane ellipse, the ellipse for the top of 
the atmosphere, and the envelope curve. 
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IX. The Shadow Forming Rays 

In all that follows, it will be assumed that the observer is on the ground, which is located 
at z = 0, and that a no-slip condition is satisfied by the mean flow. Shadow forming rays 
are those for which s z = 0 at the ground. Thus, for shadow forming rays, equation (9), 
reproduced here for convenience, requires that 


where, as indicated previously, Co is the sound speed at z = 0. Of course, only values of 
s x and s y that satisfy equation (9) are associated with rays that contribute to the shadow 
boundary. However, not all pairs (s x ,s y ) satisfying equation (9) will contribute to the 
shadow boundary simply because some of these points may be associated with rays that 
pass through a turning point between the source and the ground plane and, therefore, 
never reach the ground. Further, although the atmosphere between the source and the 
observer can filter out rays that belong to this set, and thereby remove points from the 
shadow boundary, this intervening atmosphere can never cause rays that are not members 
of this set to contribute to the shadow boundary. Therefore, the shadow boundary is 
easily determined by tracing rays associated with points (s x ,s y ) belonging to this set. 
Further, if certain members of this set are associated with rays that cannot reach the 
ground, these rays cannot contribute to the shadow boundary and there is no need to trace 
them. 

Consider a case in which the temperature decreases linearly with altitude. Then, a 
circular shadow boundary will form as shown in Figure 10. Since it is known that this 
shadow boundary is circular, it is easily detennined by tracing a single ray. However, it 
should be remembered that each point (x,y) on the shadow boundary is associated with a 
single point (s x ,s y ) in the (s x ,s y ) plane. 

If a unidirectional mean flow with a speed increasing linearly with altitude is imposed, 
some of the shadow forming rays pass through a turning point before reaching the ground 
and are, therefore, removed form the set of shadow forming rays. The rays removed from 
the shadow boundary are on the segment of the source plane ellipse along which 
dU s dc 

— • -jVT + — > 0, as discussed earlier. Eliminating these values of (s x ,s y ) from the set 
dz s t dz 

of rays traced to obtain the shadow boundary, yields the arc shown in Figure 1 1 . The 
resulting shadow boundary is shown in Figure 12. Naive geometric acoustic theory 
cannot predict the disturbance field at observer locations within the shadow region. 
Therefore, determination of the shadow boundary allows quick detennination of 
observers for which the theory cannot be applied and prevents attempts to trace rays to 
observers for which no eigenray exists, thus saving computer time. 
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Figure 10. The circular shadow boundary for a point source in a medium within which 
the temperature decreases linearly with altitude. 



Figure 1 1 . The arc of pairs (s x ,s y ) that yields the shadow boundary. 
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Shadow Boundary 
Source at x = 0, y = 0 



Figure 12. The shadow boundary for a point source in flow. 

X. Conclusions 


The characterization of the rays for a source in a stratified inhomogeneous moving 
medium and the connection of the behavior of each ray to its location in the (s x ,s y ) plane 
provides for easy identification of those rays, and the associated point (s x ,s y ), and allows 
the elimination of those rays that cannot contribute to the disturbance field at the ground 
before computer time is spent tracing them. Further, the identification of those rays that 
are members of the set of shadow forming rays, but that pass through a turning point 
before reaching the ground, makes it possible to easily and accurately determine the 
shadow boundary and thereby identify observers that cannot be connected to the source 
by any ray, thereby saving computer time. 
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